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Transition metal selenide and telluride have recently receive considerable attention due to their possible 
structural relation to ferropnictide. Pressure is often used as an efficient way to modifj^ the crystal or 
electronic structure that in many cases lead to new material states of interest. Here we search the structures 
of IrTe2 up to 150 GPa using crystal structure prediction techniques combining with ab initio calculations. 
Three new stable phases under high pressures are predicted, and their electronic structure properties, 
phonon spectra, and electron-phonon couplings are also investigated. Significant reconstructions of band 
structures and Fermi surfaces are found in these new phases. Calculated results show that while the C2/m-2 
phase has bad metal behavior and very weak electron-phonon coupling, the Pa3 and I4/mmm phases have 
relatively higher electron-phonon coupling up to ~ 1.5 and 0.7, respectively. The variable- composition 
searching have been performed, newly compounds with different stoichiometries, such as IrTe^, IrTe, and 
Ir^Te, are predicted to be thermodynamically and dynamically stable at high pressures. The pressure range 
investigated here is accessible in the diamond anvil cell experiments, thus our results might stimulate 
further experimental studies. 



In recent years, 5d transitional metal iridium ditelluride (IrTe2) has attracted much interest due to its peculiar 
structural/electronic phase transitions and emergence of superconductivity upon chemical doping or inter- 
calation^ It was reported^ that IrTe2 undergoes a structural phase transition at about 250 K from a high- 
temperature trigonal (P3ml) to a low-temperaure monoclinic (C2/m) structure, accompanied by temperature- 
dependent anomalies of electron resistivity and magnetic susceptibility. IrTe2 exhibits superconductivity when 
the structural phase transition is suppressed by chemical doping or intercalation^"^. Since the structural phase 
transition is characterized by the formation and breaking of Ir-Ir bonds along the ^-axis, the superconductivity 
can be interpreted by bond-breaking-induced superconductivity^'^^. On the other hand, it has been proposed to 
understand this phase transition down to the electron -structure level. Yang et aU suggested that the phase 
transition is a charge-orbital density wave type, and superconductivity up to 3 K sets in as soon as the density 
wave transition is suppressed by Pd intercalation and substitution. Ootsuki et aV argued that the orbitally 
induced Peierls effect plays an important role in the charge-orbital instability and superconductivity of 
Iri_^Pt;cTe2. Recent optical spectroscopic measurements and density functional theory calculations indicated 
that the transition is not driven by the density wave type instability but caused by the crystal field effect which 
further splits/separates the energy levels of Te {p^ andp-^;) and Te (pj bands^. A local bonding instability associated 
with the Te 5p states was suggested as a main factor dominating the structure and the superconductivity^'^. 

It is well known that pressure can efficiently modify the atomic and electronic structures of materials, leading to 
the formation of novel materials with unusual physical properties, and crystal structure can play a crucial role in 
the appearance of superconductivity. Recently, a pressure-induced superconductivity was reported in AuTe2 and 
its bulk superconductivity appears only in the high-pressure phase^°. Therefore, it is highly desirable to search for 
new structural phases of IrTe2 under high pressures and investigate the origin of possible superconductivity. 

In this work, we use crystal structural prediction techniques to find energetically stable/metastable structures of 
IrTe2 under high pressures up to 150 GPa. Several novel crystal structures for IrTe2 are predicted in different 
pressure regions, i.e., the Pa3 one with 4 formula units (f.u.)/cell, the C2/m-2 one with 2 f.u./cell, and the I4/mmm 
one with 1 f.u./cell. The electronic structure, phonon spectrum, and electron-phonon (EP) coupling are 
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Figure 1 | Crystal structure and equation of states of the IrTe2 phases, (a) Five types of crystallographic structures of IrTe2, and (b) enthalpy relative to 
that of the Pa3 phase vs pressure and (c) volume vs pressure curves for various structures. Insets in (c) show different type of polyhedron formed 
under pressure. 



investigated for each phase. The calculated results show that the band 
structures and Fermi surfaces (FSs) for the high-pressure phases have 
dramatically reconstructed. It is found that the Pa3 and 14/mmm 
phases have relatively strong electron -phonon coupling, which will 
be of benefit to possible superconductivity. We also have carried out 
the variable- composition search. Newly found compounds such as 
IrTes, IrTe, and Ir^Te are thermodynamically and dynamically stable 
and have unusual electronic and phonon properties. 

Results and Discussion 

Structural and Bonding analysis. First, we use crystal structural 
prediction techniques to find energetically stable structures of 
IrTe2 under high pressures. Five types of crystallographic struc- 
tures shown in Fig. 1(a) are found to be the stable/metastable 
phases with lower energies. The calculated enthalpy-pressure (AH- 
P) and volume-pressure (V-P) curves are plotted in Figs. 1(b) and (c) 
for the structures interested. From Fig. 1(b), we find that with 
increasing pressure, IrTe2 undergoes a series of structural phase 
transitions, from the C2/m to Pa3 structure at about 20 GPa, then 
to the C2/m-2 structure at about 70 GPa, and finally to the 14/mmm 



structure at about 100 GPa. The C2/m structure consists of two 
formula units, containing edge-sharing Ir-Tee octahedra. Although 
we did not get any clue in advance and "predict" the Pa3 phase from 
scratch, this structure actually has been synthesized in experiment^ \ 
and our theoretical transition pressure to it agree very well with the 
experimental value. It has a cubic lattice embodying vertex- sharing 
Ir-Te6 octahedra. An Ir ion locates at the center of the lattice and 
other Ir ions lie in the middle of the lattice edges. The C2/m-2 
structure has the same symmetry as the C2/m one, but it exhibits a 
two-dimensional (2D) characteristic in the xz plane. The Ir and Te 
ions form an unique edge-sharing Ir-Tey octahedron in the lattice. 
For the 14/mmm structure, each Ir ion is surrounded by eight Te ions, 
forming an Ir-Tes cuboid, and it consists of a stack of Ir-Tes cuboid 
arranging along the z direction. It's interesting that upon pressure, 
the coordination number of Ir ions increases continuously, forming a 
series of atomic polyhedrons of Ir-Tee Ir-Tey Ir-Teg. From 
Fig. 1(c), one sees that there is a sudden drop in volume at each 
critical pressure, indicating a first-order structural phase transition 
there. The optimized structural parameters for these five IrTe2 
phases discussed above are hsted in Table I. 
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Table I | Opti mized structural parameters of five structures for lrTe2 calculated with PBE and DFT-D2 methods. The parameters of C2/ m and 
Cl/m-2 shown here correspond to the primitive cell. Note that the volumes are unified in one formula unit 
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As found in many other system, such as CO^^ and C02^^, some of 
the IrTe2 phases present a laminar crystal structure, we did additional 
calculations including van der Waals (vdW) correction proposed by 
Grimme (DFT-D2)^^. The obtained structural parameters are listed 
in Table I. The message is that the volume of unit cell for several 



structures shrinks due to the vdW corrections. The C2/m and P3ml 
structures under ambient condition have the largest percentage of 
contraction (about 5.3%) and we find this is mainly from the decreas- 
ing of the distance between the adjacent IrTe2 layers. However, the 
magnitude of contraction decreases with increasing pressure, reflect- 




Energy (eV) Energy (eV) 



Figure 2 | Electronic structures of the novel IrTe2 phases. Total and partial density of states, fat-band structures with Ir-d and Te-p characteristics, and 
Fermi surfaces for C2/m at ambient pressure (a) - (c), Pa3 at 40 GPa (d) - (f), Cllm-2 at 80 GPa (g) - (i), and 14/mmm at 120 GPa (j) - (1). 
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Figure 3 | Calculated phonon dispersions, projected phonon densities of states (PDOS) and Eliashberg spectm (from left to right) (a) C2/mat 

ambient pressure, (b) Pa3 at 40 GPa, (c) Cllm-2 at 80 GPa, and (d) lAlmmm phases at 120 GPa. 



ing that the vdW contribution under high pressure becomes small. 
Based on the electron locaUzation function (ELF), as showned in the 
Supporting Informtion, we find there are small overlap between the 
electrons on the Te atoms in two adjacent layers, both in P3ml and 
C2/m structure, this somehow gives the reason of why the binding 
energy between the layers are quite large. Without vdW corrections, 
we obtain the binding energies are 219 meV for P3ml and 214 meV 
for C2/m structure. While with vdW corrections, we obtain 847 meV 
and 843 meV, respectively. 

The calculated total and partial density of states (DOS), band 
structure with Ir-d and Te-p characteristics, and Fermi surfaces are 
shown in Fig. 2. For the DOS of the Cllm structure, the Ttp^ orbital 
exhibits a distinct peak located at the Fermi level {Ep), indicating a 
small van Hove singularity (vHs) behavior. A recent comprehensive 
study of the electronic states has revealed the vHs origin of the phase 
transition between high temperature P3ml phase and low temper- 
ature Cllm phase in IrTea^^ According to the orbital characteristics 



of the fat-band, two hybridized Ir 5d and Te 5p bands cross the Fermi 
level, forming two Fermi surfaces. Compared with the P3ml struc- 
ture, the interlayer interaction of the C2/m structure and hence its 
band dispersion along the c axis (F-A direction) are suppressed, so 
that the 2D character of the FSs becomes enhanced. The reconfigura- 
tion in the topology of FSs along the direction appears to be 
consistent with what has been observed in angle-resolved photoemis- 
sion spectroscopy measurements^^"^^. 

The Pa3 structure exhibits quite unique electronic structure and 
Fermi surfaces. There are four bands crossing the Fermi level, which 
construct four compUcated 3D FSs. A van-Hove-like peak Ues below 
the Fermi level 0.1 eV, yielding a large DOS at the Fermi level (N(Ep)) 
~ 7 states/eV per unit cell, which is contributed mainly by the Ir 5d 
planar orbitals and the Te p^ orbital. A pseudo gap appears around E 
= —1.9 eV, where a significant Fermi cone originates mainly from 
the Tep orbitals around the T point in the band structure. The Fermi 
level of the C2/m-2 structure lies at the valley bottom of the DOS 
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curves and N(Ep) ~ 1 states/eV per unit cell, exhibiting a bad metallic 
behavior. Two bands crossing the Fermi level form the FSs shown in 
Fig. 2(h). For the 14/mmm structure, the FSs arise from the four 
bands crossing Ep. Around the T point there is a claw-like Fermi 
pocket with four-fold symmetry, and the Fermi pockets at zone 
corners are shaped as irregular cylinders. At the same time, there is 
a reduction of DOS at Bp for 14/mmm, N(Ep) = 1.33 states/eV, 
compared with those in the low-pressure phases. The reconstructed 
band structures and FSs may greatly affect phonon properties and 
electron-phonon coupling. 

Vibrational Analysis. The calculated phonon dispersions, projected 
phonon densities of states (PDOS), and Eliashberg spectral function 
oc^F{w) are shown from left to right in Fig. 3. In the left panel of Fig. 3, 
the distribution of the EP coupling for each phonon mode, Iqy, is 
described by red circles, the radius of the circles being proportional to 
the magnitude of each >^^^. For the C2/m structure, as shown in 
Fig. 3(a), the two gaps near 120 cm"^ and 150 cm"^ in the phonon 
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Figure 4 | Calculated superconducting properties of IrTe2 phases as 
functions of pressure, (a) Calculated electron-phonon coupling and 
(b) superconducting transition temperature. 



spectrum divide the spectrum into three regions: the acoustic 
branches, middle-frequency and high-frequency optical branches. 
The acoustic modes extend to about 100 cm~\ The modes at the 
high frequency region are associated with the vibrations of Ir atoms 
beating against Te atoms. Whereas the modes at the middle 
frequency region are associate with the vibrations of Te-Te bond 
stretching. The obtained a^F{co) curve exhibits two main humps 
below 120 cm~\ which originate from the acoustic modes. 

For the Pa3 structure at 40 GPa [see Fig. 3(b)], phonon frequen- 
cies extend up to ~ 280 cm" ^ and the point group at the T point is T/^, 
which can be decomposed as F = 0 £^ 0 3T^ 0 2A^ 0 lE^ 0 
6Tu with Eg and E^ being double degenerate, and Tg and T„ triple 
degenerate. There are two tiny gaps in the PDOS, respectively, near 
125 cm"^ and 210 cm"\ The acoustic phonon modes exhibit a sig- 
nificant phonon softening along all the symmetric directions. 




Figure 5 | Formation enthalpies of Ir^Tei_^per atom with respect to Ir 
and Te for different Ir-Te phases. The abscissa x is the fraction of Ir in the 
structures. The blue circles indicate the most stable structures at each 
stoichiometry. 
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Table 3 Optimized structural parameters of IrTea, IrTe and IraTe, calculated 
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Furthermore, we notice that the optical branches at the lower fre- 
quency side also become soften and merge with the acoustic 
branches. In the DFPT calculations, the Eliashberg spectral function 

depends directly on the EP matrix element, k;^km(k+q) ' which is 

defined as the variation of the self- consistent crystal potential. The 
low-frequency spectral function curve (X^F(w) has two peaks near 
50 cm"^ and 100 cm"\ which are contributed by the acoustic and 
optical T„ modes. The latter originate from the out-plane breathing 
and rocking vibrations of the Ir and Te atoms. The integration of 
X{cd) up to 100 cm"^ is near 0.7, almost 80% of the total EP coupling. 
In the Pa3 phase, the softening of low-frequency modes and 
enhancement of the Eliashberg spectral function can significantly 
increase the EP coupling. 

For the C2/m-2 structure at 80 GPa [see Fig. 3(c)] , the point group 
at the r point is Q, same as the C2/m phase, which can be decom- 
posed as r = 9Au 0 9Ag. Both acoustic and optical branches are 
hardened and the obtained EP coupling strength is very weak. For the 
14/mmm at 120 GPa [see Fig. 3(d)], the phonon frequencies expand 
up to more than 300 cm"^ and the point group at the T point is 04^. 
The r modes can be decomposed as F = 2A2„ 0 0 £^ 0 Aj^, 
where and Eg modes are double degenerate. From the PDOS 
curves, one sees that the low-frequency branches below 125 cm"^ 
are mainly associated with the vibrations of the Ir atoms, whereas 
the high-frequency branches above 200 cm"^ originate mainly from 
the vibrations of the Te atoms. The (X^F(w) curve exhibits a distinct 
peak around w = 100 cm"\ and the low- frequency phonons 
account for more than 60% of the total EP coupling. The calculated 
phonon vibration frequencies for each IrTe2 compound at the T 
point are Usted in Table 2. 

Using the Allen-Dynes formula^^'^^ 



T.-^exp 



1.04(l + >^) 



2-Ai*(l + 0.622) 
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where co/„ is the logarithmically averaged frequencies, and taking 
the Coulomb parameter as = 0.15, we evaluate the pressure- 
dependent Tc. The obtained X and as a function of pressure are 
plotted in Figs. 4 (a) and (b), respectively. It is found that the calcu- 
lated value is equal to /I ~ 0.4 for the C2/m phase and X ~ 0.25 for the 
C2/m-2 phase. Such weak EP couplings, together with bad metallic 
behavior in the C2lm-2 phase, result in almost vanishing T^. Very 
interestingly, the predicted value of k for the Pa3 phase is enhanced 
up to ~ 1.5, yielding high-pressure induced superconductivity with 
Tc ~ 6 K. In the 14/mmm phase, the maximal X is equal to 0.7 and it 
decreases with further increasing pressures. 



Variable-composition Search. We also investigated the phase 
stabilities of Ir-Te systems calculating the formation enthalpy of 
various Ir^cTei-;^ compounds, as shown in Fig. 5. The formation 
enthalpy £yof Ir^^Tei-^^ was calculated as 

Ef{Ir,Tei_,) =E{lr,Tei_,) -xE{Ir) - (1 -x)E{Te) (2) 

In the entire region of explored pressures, IrTe2 stays as a stable 
compound. The Ir-rich phase emerges at high pressure. Convex 
hull diagram in Fig. 5 features newly predicted thermodynamically 
stable compounds: IrTca at 20 GPa [Trigonal (space group R3c)]; 
IrTe at 70 GPa [Monoclinic (space group C2/c)] and 120 GPa 
[Orthorhombic (space group Cmcm)]; and Ir3Te at 120 GPa 
[Tetragonal (space group 141 mem)]. Electronic densities of states 
show that all of these compounds are good metals. The lattice 
parameters including the vdW (DFT-D2) correction are shown in 
Table 3. In IrTca and two phases of IrTe, six Te atoms bonding with 
one Ir atom construct a Ir-Tce octahedra. In IraTe, the fundamental 
unit is consisted by Te-topped cube with Ir atoms locating at the face 
centers (Fig.Sl in supplementary information). For all the structures, 
we calculated phonons and EP couplings, and found all of them to be 
dynamically stable (Fig.S2). The EP couplings given by phonon 
calculations are from 0.3 to 0.5. 

In summary, with crystal structural prediction techniques, it is 
found that, with increasing pressure, IrTe2 may undergo a series of 
structural phase transitions: C2 / m^Pa3 -^C2/ m-2^14 jmmm, and 
the latter three phases under high pressures are newly discovered. For 
various phases of IrTe2, the electronic band structures, phonon spec- 
tra, and electron-phonon couplings have been studied. It is shown 
that there are relatively higher EP coupling constant in the high- 
pressure Pa3 and 14/mmm. We hope such a pressure induced 
superconductivity in IrTe2 will stimulate new experiments on this 
material. The variable-composition search gives unexpected chem- 
ical compounds at high pressure with good metal feature. These 
results suggest that pressure is an effective way to overcome the 
kinetic barrier to formation in the synthesis of Ir-Te binary system. 

Methods 

Two crystal structure prediction methods, crystal structure analysis by particle swarm 
optimization ( CALYPSO) ^""^^ and ab initio random structure searching (AIRSS)^^'^* 
are employed. Both methods give the same best candidates independently for given 
chemical compositions and external conditions. The electronic structure calculations 
with high accuracy for the stable structures are performed using the full-potential 
linearized augmented plane wave (FP-LAPW) method implemented in the WIEN2K 
code^^. The generalized gradient approximation (GGA)^^ is applied to the exchange- 
correlation potential calculation. The muffin tin radii are chosen to be 2.5 a.u. for Ir 
and 2.33 a.u. for Te. The plane-wave cutoff is defined by RKj^^x = 7.0, where R is the 
minimum LAPW sphere radius and K^ax is the plane-wave vector cutoff. The pho- 
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non calculations are carried out by using a density functional perturbation theory 
(DFPT) approach through the Quantum-ESPRESSO program^^. The cutoffs are 
chosen as 30 Ry for the wave functions and 300 Ry for the charge density. The 
electronic integration is performed over al6X16X16 /c-point mesh. Dynamical 
matrices and the EP couplings are calculated on a 4 X 4 X 4 ^-point grid. A dense 24 
X 24 X 24 grid is used for evaluating an accurate EP interaction matrix. Due to the 
large atomic numbers of Ir and Te, the spin-orbit coupling is included in our calcu- 
lations. 
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